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_ ., ■ Abstract 

A Neutrino Factory producing an intense beam composed of u e (js e ) and 
>*n ; v^Vp) from muon decays has been shown to have the greatest sensitivity 

to the two currently unmeasured neutrino mixing parameters, #13 and 8c p. 
Using the 'wrong-sign muon' signal to measure v e — > v^(V e — > V^) oscilla- 

^ ! tions in a 50 ktonne Magnetised Iron Neutrino Detector (MIND) sensitivity 

to 5c p could be maintained down to small values of #13. However, the de- 
tector efficiencies used in previous studies were calculated assuming perfect 

qq ■ pattern recognition. In this paper, MIND is re- assessed taking into account, 

for the first time, a realistic pattern recognition for the muon candidate. Re- 
optimisation of the analysis utilises a combination of methods, including a 

^ \ multivariate analysis similar to the one used in MINOS, to maintain high 

O efficiency while suppressing backgrounds, ensuring that the signal selection 

efficiency and the background levels are comparable or better than the ones 
■ in previous analyses. 
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1. Introduction 

The concept of a neutrino beam from the decay of muons in a storage ring 
was first proposed in 1980 [1]. More recently, such a facility was explored 
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as a preliminary stage towards a muon collider and was renamed "Neutrino 
Factory". Its physics potential was originally described by Geer jij. The 
great advantage of a Neutrino Factory over conventional neutrino beams 
from pion decay is that the decay of muons is very well described by the 
Standard Model and so the beam flux is easily calculable. Therefore, it is 
possible to perform high precision neutrino oscillation experiments at a high 
flux Neutrino Factory. Another significant feature of a Neutrino Factory is 
that one can accelerate muons of both signs into a storage ring, thereby 
enabling study of both neutrino and anti-neutrino oscillations with equal 
flux, vastly improving sensitivity to CP violation in the neutrino sector. For 
a more recent review see 0] . 

Early papers on the physics outcomes of a Neutrino Factory concentrated 
on the sub-dominant v e — > oscillation 0] in which a muon of opposite 
charge to that stored in the facility storage ring (wrong-sign muon) would 
be produced in a far detector by the charge current (CC) interactions of 
the oscillated v^. The first analysis of the capabilities of a large magnetised 
iron detector to detect the wrong-sign muon signature was discussed in j^j 
(termed the Golden Channel), where it was demonstrated that this combi- 
nation was capable of the extraction of the remaining unknown parameters 
in the neutrino sector, the third mixingangle #13 of the Pontecorvo-Maki- 
Nakagawa-Sakata (PMNS) matrix 0, 0,0] and the CP violating phase Sep- 

The Magnetised Iron Neutrino Detector (MIND) is a large scale iron and 
scintillator sampling calorimeter. As a result of the studies mentioned above 
it is considered the baseline detector for a Neutrino Factory (NF) sto ring 
muons in the energy range 20-50 GeV 0] . Under the remit of EUROnu fief 
and the International Design Study for a Neutrino Factory ll[ all aspects of 
possible future neutrino beam facilities including accelerator, detectors and 
physics must be studied and compared to select the best option to determine 
the remaining oscillation parameters. 

Previous studies of MIND focused on the topology and kinematics of neu- 
trino events in the detector, assuming perfect pattern recognition. By smear- 
ing the kinematic variables of the participant muon and hadronic shower it 
was demonstrated that using a combination of cuts on the relative length of 
the two longest particles in the event and the momentum and isolation of this 
candidate, high si gna l identification efficiency and background suppression 



could be achieved [12], |13| . However, a full study without such assumptions 
is necessary to fully characterise the detector response. 

While MIND is essentially a large scale version of the MINOS detec- 
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tor jl4{, the nature of the NF beam - containing 50% u e and 50% in the 
case of stored fi + - means that the optimisation of the analysis is some- 
what different. Incorrect charge assignment (charge misidentification) of 
non-oscillated CC interactions present a significant possible background 
in this beam configuation, in addition to backgrounds from meson decays in 
the hadronic shower and misidentification of Neutral Current (NC) and v e 
CC events. 

This current study re-visits the problem by taking an un-biased look 
at the visible part of a large sample of neutrino interactions - generated 
using the same GEANT3 15] simulation as in the above mentioned studies 
with a uniform distribution in neutrino energy - and developing pattern 
recognition algorithms (first presented in [l6j]) - described in Sec. [3] - to 
extract a candidate muon for fitting using a Kalman filter. Successful fits 
are then subject to offline analyses - described in Sec. H]- to determine the 
validity of those wrong sign candidates. Analysis results are presented in 
Sec El 



2. MIND parameterisation and expected event yields 

For the purpose of the described study, MIND is a cuboidal detector 
of 14 m x 14 m cross-section and 40 m length, segmented as 4 cm of iron 
and 1 cm of plastic scintillator for a total mass of ~51.0 ktonnes. A dipole 
magnetic field of mean induction 1 T in the transverse plane provides the 
field neccessary for charge and momentum measurements. 

In the first part of the analysis, event vertices were generated centred 
in the detector plane at 1.5 m from the front of the detector in the beam 
direction (z) in order to study the nature of the backgrounds without detector 
edge effects. Sec. 14.31 discusses the expected fiducial effects when a more 
realistic randomly generated vertex is considered. 

At a MIND placed 4000 km from the neutrino source and assuming the 
current best global fit oscillation parameters: 9u = 33.5°, #13 = 5.7°, 6*23 = 
45°, Am 2 21 = 7.65 xl0~ 5 eV 2 , Am 2 32 = 2.40 x 10~ 3 eV 2 [l7|, setting 5 C p = 45° 
and calculating matter effects using the PREM model [18| , the expected total 
number of interactions due to 10 21 /i + decays at 25 GeV energy would be of 
order those shown in table [TJ 

Thus in order to successfully extract oscillation parameters from the 
golden channel, potential backgrounds from non-signal interactions must be 
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v e CC 


77^ + v e NC 


(Signal) 


1.22 x 10 5 


3.34 x 10 5 


1.48 x 10 5 


5.56 x 10 3 



Table 1: Expected absolute number of interactions in a 51 ktonne MIND at a distance of 
4000 km from a NF storage ring with 25 GeV muons. 

suppressed to at most the 10 -3 level in absolute terms. Moreover, the ex- 
istence of possible degenerate solutions due to uncertainty in the measured 



parameters and due to the nature of the oscillation probability (see |l9l. l20l|) 
means that spectral information is required to reliably determine 5c p. This 
additional requirement dictates that backgrounds must be suppressed to be- 
low 10~ 3 in each energy bin while maintaining an efficiency threshold below 
5 GeV so that information on the rise of the first oscillation maximum is 
available. 

3. Reconstruction tools 

The reconstruction package was used to analyse a large data set com- 
prised of Deep Inelastic Scattering (DIS) neutrino interactions of and u e 



generated by the LEPT061 [21J package and tracked through the GEANT3 
simulation of MIND. Considering CC interactions of and v e with a dedi- 
cated study of events containing wrong sign muons from meson decay in 
CC and NC interactions the main expected backgrounds were studied. 

Each event considered comprised all three dimensional points with their 
associated energy deposit, which were recorded in the scintillator sections of 
the MIND simulation, with the x,y position of these hits smeared according 
to a a — 1 cm Gaussian before analysis began. 

3.1. Muon candidate extraction 

After ordering the hits from smallest to greatest z position in the detector 
the first step of the reconstruction was to extract a candidate muon from the 
event. Two methods were employed to perform this task depending on the 
event topology: a Kalman filter incremental fit was used to extract candidates 
from those events with one particle clearly longer than the others (described 
in Sec. 13 . 1 . 1 [) . while a Cellular Automaton method was used in those events 
not viable for reconstruction through the first method (see Sec. I3.1.2p . The 
criterion for the first category was that the five planes with activity furthest 
downstream should contain no more than one hit per plane. 



4 




Figure 1: Muon candidate purity as a function of true muon momentum for i) the Kalman 
filter extraction method and ii) the Cellular automaton method. 



3.1.1. Kalman Filter candidate extraction 

Using the Kalman filter algorithm provided by RecPack 22] it is possible 
to propagate the track parameters back through the planes using a helix 
model, which takes into account multiple scattering and energy loss. Since, 
in general, a muon will act as a Minimum Ionising Particle (MIP) and will 
travel further in the detector than the hadronic particles, those hits furthest 
downstream can be assumed to be muon hits and used as a seed for the 
Kalman filter. The seed state is then propagated back to each plane with 
multiple hits and the matching x 2 to each of the hits is computed. Hits 
with matching \ 2 below 20 are considered and in each plane the one with 
the best matching among these is added to the trajectory and filtered (the 
track parameters are updated with the new information). All accepted hits 
constitute the candidate muon and are presented for fitting (Sec. 13. 2p . with 
the remaining hits being considered as hadronic activity. Fig. [T]-(left) shows 
the fraction of true muon hits in the candidate when using this method. 

3.1.2. Cellular Automaton candidate extraction 

Events with high Q 2 transfer or low neutrino energy can tend to be re- 
jected by the first method, since in general the muon will not escape the 
region of hadronic activity. In order to recover these events a second method 
is employed. The Cellular Automaton method (based on the method de- 



scribed in [23|) uses a neighbourhood function to first rank all the hits and 
then form all viable combinations into possible trajectories. 

A 'neighbour' is defined as a hit in an adjacent plane within a pre-defined 
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transverse distance of the projection into that plane of the straight line con- 
necting hits in the previous two planes. Starting from the plane with lowest 
z position, hits are given a rank one higher than their neighbour in the pre- 
vious plane should they have one. Trajectories are then formed from every 
possible combination of one hit per plane starting with those of highest rank 
using the neighbourhood function with a stricter condition. 

Those trajectories formed using this method are then subject to a number 
of tests to determine which is most likely to be a muon. After having a basic 
helix fit performed and being assessed according to their length, trajectories 
are rejected for being short, having high \ 2 fit or high relative curvature 
error (described in Sec. 14.11) . The candidate muon is then selected as the 
longest remaining trajectory with the lowest y 2 . All other hits in the event 
are considered to be from hadronic activity. Fig. [1]- (right) shows the purity 
of the candidate when using this method. 

3.2. Candidate fitting 

All candidates successfully extracted from their event that have greater 
than six hits are presented to the fitter as a candidate muon. The same 
Kalman filter algorithm is used here as in Sec. 13.1.11 Fitting the candidate 
iteratively improves seeding and thus using a more constricted \ 2 condition 
than in the pattern recognition, the maximum number of successful, reliable 
fits were achieved. 

With the trajectory hits ordered in increasing z position, a least squares 
quartic fit was performed on the section outside the planes where there was 
hadronic activity. This fit was used to estimate the slopes in x and y and the 
momentum of the candidate, to be used as a seed for the Kalman filter helix 
fit in the forward direction. The matching x 2 was once again checked at each 
hit. Hits with greater than the pre-determined maximum (20) were ignored. 
In addition, the filtering process only allows a pre-determined maximum 
number of hits (5) to be ignored. Should this number be reached, the filtering 
process is aborted and the smoother uses only those hits up to this point in 
the candidate. This method efficiently rejects hits beyond any large angle 
scatter which could cause charge misidentification. Successful fits were re- 
seeded with the state vector at the first fitted hit and a scalar multiple (5) 
of the corresponding covariance matrix (taking only the diagonal elements) 
and then refitted. 

Failed fits and those with less than 50% of their hits fitted are then fitted 
again in the backwards direction using the seed from the pattern recognition. 
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Two iterations are once again performed, with successful fits being accepted 
and those which are unsuccessful reverting to the result of the original fit. 

The result of a fit being the track parameters at the projection to the 
event true vertex z position (3- momenta, position and charge). 

3.3. Hadronic reconstruction 

The hadronic activity must be used to reconstruct the energy of the 
hadronic shower in order to ultimately reconstruct the energy of the inter- 
acting neutrino. In the absence of a well developed algorithm to perform 
this task, the current study assumes reconstruction of the hadronic energy 
Ehad with a resolution SE had equal to that recorded by the MINOS CalDet 



testbeam [14|, |24 



^ = 4^1 ©0.03. (1) 

Ehad V Ehad 

It was demonstrated in [5| that a cut based on the isolation of the muon 
candidate from the hadronic shower was a powerful handle for the rejection 
of hadronic backgrounds. This isolation was measured via the Qt variable: 

Qt = P, sin 2 0, (2) 

where is the muon momentum and is the angle between the muon 
and the resultant hadronic vector. This requires the reconstruction of the 



direction vector of the shower. The Monolith test-beam 25] measured an 
hadronic angular resolution described by: 

10 - 4 w - 1 
56 h ad = -j== © p— (3) 

V t^had ^had 

for a similar detector. This parameterization was used to smear the hadrom 
shower direction vector, which in combination with the reconstructed muon 
momentum and direction (see Sec. 13.21) were use to compute the Qt variable 
defined above. 



4. Analysis tools and cuts 

As mentioned in Sec. [1] there are four main possible sources of background 
to the wrong sign muon search: incorrect charge assignment and high energy 
wrong sign muons from meson decays in F M CC events, and NC and v e CC 
events wrongly identified as CC. In order to reduce these backgrounds 
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while maintaining good efficiency a number of offline cuts were employed. 
They can be organised in four categories: i) muon candidate quality cuts, ii) 
CC selection cuts, iii) fiducial cuts and iv) kinematic cuts. 

4-1. Muon candidate quality cuts 

These cuts are related to the quality of the candidate track fit and the 
determination of its curvature. Two observables are considered: the x 2 prob- 
ability of the Kalman filter fit and the relative error of the determined curva- 
ture (■dzf )■ The first helps in rejecting high angle scatters or muon candidates 
with a large contamination from hadronic hits. The second variable is related 
with the probability of misidentifying the charge, and shows significant sepa- 
ration for correct and incorrect charge assignments as shown in Fig. [2]-(left). 

It is possible to reject a large portion of possible backgrounds using se- 
quential cuts on these two variables: 



O-q/p 



q/p 



< 0.7 and xlrob > 0.9999, (4) 



(where Xprob * s the probability as calculated in the TMath class of the 
ROOT framework 26]). However, a slightly better rejection is found when 
the relative error cut is substituted by a cut on the log likelihood ratio of 
-H^- for signal and background (see Fig. [2]-(right)): 

C q/P > 2 and xlrob > 0.9999. (5) 



4-2. CC selection cuts 

The discrimination between CC and NC interactions relies on three 
easily available or calculable parameters, which are those of the extracted 
muon candidate. Due to the similarity of MIND and MINOS the parameters 
employed in the MINOS analysis 23] were used. Using a high statistic data 
set with knowledge of the true nature of each event, distributions of these 
three parameters for both NC and CC events were formed into PDFs (or 
likelihoods). 

The first parameter was the length of the candidate in terms of the number 
of hits which form it (Ihu)- This variable takes advantage of the nature of 
the muon as a penetrating particle and shows clear separation between 
CC and NC events (see Fig. H (top-left)). 
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til background from v„ GC 



5 10 

Log likelihood of fractional error 



Figure 2: i) SiE. likelihood for signal and background from v„ CC, and ii) q [ v log 

q ' p '' 
likelihood ratio (C q / V ). 



The second parameter is the fraction of the total visible energy in the 
event which is in the candidate (Z/ rac ). This parameter is not useful for 
all events due to the high probability for both NC and CC events to have 
a fraction very close to or equal to one. Thus, events that fall into this 
category, low Q 2 CC events or single pion production NC predominantly, 
are excluded from this distribution and do not use this parameter in their 
analysis. Here, while NC events demonstrate the full spectrum of possible 
values, signal events tend to be more concentrated at high fractions (see 
Fig. [3]- (top- right)). However, high Q 2 CC events will tend to exhibit NC like 
behaviour. 

While the third parameter used by MINOS is the mean energy deposited 
per plane for the candidate, the current simulation setup of MIND does not 
exhibit sufficient separation in this parameter for effective analysis. Thus, in 
place of this parameter the variance of the deposit is used (l var ), shown in 
Fig. El- (bottom). 

The likelihood ratio for each of the three observables was computed and 
combined in three main log likelihood discriminators described in Eqs. [6]to[8j 

c 2 



( ]CC v ice v ]CC \ 

, / L hit A L frac A L var \ /„n, 

° g 1 NC X 1 NC x Wc ^ 

\ l hit A l frac A L var J 
/ JCC x JCC \ 

1 f hp x L hu \ (7) 

iU & 1 ]NC v ]NC ) \' 1 

\ l 2D x L hit J 
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Figure 3: PDFs of the three parameters used for NC/CC likelihood separation, i) Number 
of hits in candidate Ihu, H) Fraction of visible energy in candidate lf rac an d Hi) Variance 
of energy deposit in candidate l var . 



10 




— v (l CC events 
---- NC events 




log likelihood for 2 parameter mult. rat. ('single track') 



Figure 4: Log likelihood discriminator distributions for energy fraction < 0.999, i) C\, 
defined in Eq. ii) £2, defined in Eq. [7| and Hi) for energy fraction > 0.999, £3, defined 
in Eq. 

(j CC frac=l .Q C frac=l \ 

l hit 2 W \ (o\ 

lNC frac=l 1NC ; rac=1 I \° ) 

l hit A l var J 

C\ is formed from the multiplication of the likelihoods mentioned above while 
£2 is formed by the multiplication of the lut likelihood and a 2 dimensional 
likelihood of the variance and energy fraction (1 2 d = lhit '■ Ifrac)- £1 or £-2 
are used when the energy fraction is less than 0.999 and £3 otherwise. Dis- 
tributions of these discriminators for samples of NC and CC events are 
shown in Fig. HI 

4-3. Fiducial cuts 

Events originating near the edges of the detector can leave the sensitive 
volume. This will result not only in the loss of event energy and thus wors- 
ened energy resolution but, due to the shortening of the event, can cause 
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a misidentification of the charge of a candidate. While the shortening of 
the event has the potential to reduce backgrounds from NC and u e CC as 
there should be less viable candidates, viable signal can also be lost, with 
a corresponding increase in charge misidentification background. Therefore, 
it is recomendable to apply a fiducial volume cut so that these pathologies 
are minimised. Specifically, events are rejected should their candidate have 
both its first hit within 50 cm of the sides or back of the detector and its 
last within 10 cm. In sec. 15.1.31 the edge effects and their suppression are 
presented using z/ M CC events as a model since they should affect little or no 
increase on NC and v e CC backgrounds and any small variation should be 
of the same spectral form as those seen in CC events. 

4-4- Kinematic cuts 

Considering the remaining signal and background after applying all cuts 
described above, the Q t (see Sec. 13. 3f) and muon candidate momentum (P M ) 
distributions are those shown in Fig. |5j A clear separation between signal 
and background events is observed. In particular, background events are 
concentrated at very low Qt, while the signal exhibits much larger Q t values. 
In order not to reduce the efficiency at low neutrino energy, cuts on these 
two variables are only applyed for reconstructed neutrino energy (E u ) above 
7 GeV. The applyed cuts are those of Eq. [9j 

P M ^ 0.2 ■ E v and Q t > 0.25 GeV/c for E v > 7 GeV. (9) 

4-5. Summary of analysis cuts 

As will be dicussed in the next section the most succesful set of cuts is 
given below in table EJ 

5. Analysis Results 

Using a large data set and the analyses described above the efficiency and 
rejection power of MIND has been studied. 

5.1. Up charge current interactions 

The background from interactions can be separated in two different 
contributions: i) fake wrong-sign muons from charge misidentification of the 
primary muon (mainly) and from pion to muon confusion, and ii) true wrong- 
sign muons from the decay of hadrons. 
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Figure 5: Distributions of kinematic variables: Reconstructed muon momentum (left) and 
Q t variable in GeV/c units (right) versus reconstructed neutrino energy, for (top— > bottom) 
signal, overlinev^ CC backgrounds, NC and v e CC. 



Cut type 


Cut value 


Fiducial 


if \rj irst — rf et \ ^ 50 cm, r' as * — rf et \ > 10 cm 
for Ti — x,y,z 


Track quality 


C q/p > 2.0 and X 2 P rob > 0-9999 


Ufj, CC selection 


Ci > for Z/roc < 0.999 
£ 3 > for l frac > 0.999 


Kinematic 


> 0.2 ■ E„ and Q t > 0.25 for E v > 7 GeV 



Tabic 2: Summary of analysis cuts. 

5.1.1. Incorrect charge assignment 

The charge misidentification background was studied using interac- 
tions where events containing hadronic decays to /i~ were excluded to be 
considered separately (Sec. 15.1. 2p . An event is considered background if a 
candidate is successfully extracted and fitted with charge opposite to that 
of the true primary muon. Background events are mainly due to incorrect 
charge assignment to the true primary muon (due to multiple scattering or 
impurity of the candidate), but have a small contribution from penetrating 
hadrons (mainly pions) which are identified as muon candidates when the 
true primary muon has low momentum and is not correctly identified. 

As shown in Fig. [6]-(left) this background can be efficiently suppressed by 
cutting on the track quality variables, described in Sec. 14.11 Further rejection 
is obtained by applying CC selection cuts (see Fig. El-(right)). 

5.1.2. Wrong sign muons from hadron decays 

The production and decay of negatively charged mesons in the hadronic 
part of a DIS interaction has high probability to produce a fi~. Particu- 
lary mesons containing charm will decay promptly and produce high energy 
muons which can be selected as primary muon candidates when the true 
primary muon is not correctly identified (in general when it has low mo- 
mentum). Suppression of this background is particulary important due to 
the high level of uncertainty on the value of the charm production cross sec- 
tion [28(. Track quality and CC selection cuts are effective in reducing 
this type of background for low neutrino energies, but the suppresion of high 
energy background requires further cuts, which are described below. 
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Figure 6: Charge mis-assignment and hadron to muon mis-identification background i) 
with track quality cuts only and ii) including CC selection cuts. 
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Figure 7: Wrong sign muon decay background with track quality and CC selection cuts. 
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5.1.3. Inclusive background with fiducial and kinematic cuts 

Considering a set of z7 M CC events generated randomly throughout the 
detector volume the inclusive background from this type of interaction has 
been studied. An additional cut on those events with a candidate failing the 
fiducial volume cut, defined in Sec. 14.31 is used to suppress background caused 
by edge effects. As can be seen in Fig. IB]- (left), the additional background 
introduced by edge effects is almost compensated by the fiducial cut, leading 
to a inclusive I7 M background similar to the addition of the ones shown in 
Figs. [6]- (right) and [7] 

The high energy background of Fig. (B]-(left) is mainly due to very hard 
muons from the decay of charm mesons. Fortunately these muons, due to 
their decay origin, tend to be imbedded in the hadron shower unlike a true 
primary muon. Thus the Qt variable should be very effective in rejecting this 
kind of event. As can be seen in Fig. IB]- (right) the kinematic cuts afford a 
sizeable suppression, particularly at higher neutrino energy. 




Figure 8: Expected background from v ^ CC interactions when events are randomly gener- 
ated in the entire detector: i) after track quality, CC selection and fiducial cuts, and 
ii) including kinematic cuts. 



5.2. Neutral current interactions 

Neutral current interactions should be of the same nature for all species. 
Background events will tend to originate from penetrating pions or muons 
from the decay of hadrons. Moreover, since there will always be missing en- 
ergy in the event, those events successfully fitted will tend to be reconstructed 
at lower energy than the true neutrino energy. As such and due to the large 
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Figure 9: Expected background from NC interactions, i) with track quality cuts (Eq.^ 
only, ii) including CC selection cuts, Hi) substituting track quality cuts of Eq. ^ by 
those of Eq. [5| and iv) including kinematic cuts. 



amount of NC events expected in the detector, this background must be sup- 
pressed efficiently. Fig. M shows the evolucion of the NC background when 
different cuts are included. 

5.3. u e charge current interactions 

The interactions of v e present in the beam can also produce some back- 
ground to the signal. While the electron itself will be stopped quickly and will 
shower far more than a muon, penetrating pions or decay muons originating 
in the hadronic shower can be mistaken for primary muons. 

After application of track quality and CC selection cuts it can be seen 
in Fig. [10]- (left) that the v e CC background can be reduced to a similar level 
as the NC background when the same cuts are applied. Moreover, excellent 
rejection for high energy neutrinos is obtained by applying kinematic cuts, 
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Figure 10: v e CC background to golden channel signal, i) with track quality and CC 
selection cuts, and ii) including kinematic cuts. 



as shown in Fig. [TOUright). This is because the candidate muon in v e events 
tends to have lower momentum than in NC events, as shown in Fig. [5] . 

5.4- Summary 





Total 


0-5 GeV 


5-10 GeV 


10 - 30 GeV 


v„ CC 
v e CC 
+ u e NC 
(signal) 


5.5 x 10" 4 
7.8 x KT 6 
3.8 x 1(T 5 
0.64 


6.6 x 10" 4 
2.6 x i0~ 5 
6.8 x 10~ 5 
0.25 


8.2 x KT 4 
5.5 x 10~ 6 

4.3 x lO" 5 
0.66 


4.6 x 10" 4 
5.5 x 10~ 6 
3.1 x KT 5 
0.69 



Table 3: Summary of expected fractional signal and background with true neutrino energy. 

Considering all types of events mentioned above and applying the most 
successful analysis chain described in table |2] the resulting signal efficiency 
and fractional backgrounds are those summarized in table [31 The evolution 
of the backgrounds for the different cuts is in Figs. [H [9] and [TOj while similar 
plots for the signal efficiency are those of Fig. [TTJ While it is obvious that the 
effect of the kinematic cuts below 7 GeV true neutrino energy is small due to 
their application only to events reconstructed with energy greater than this 
value, it is important to remark that fiducial cuts do not affect the efficiency 
at low energies either. 
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Figure 11: Expected signal identification efficiency: i) after track quality and CC selec- 
tion cuts, ii) including fiducial cuts, and Hi) including kinematic cuts. 



Another important question is that of the relation between the true 
and the reconstructed neutrino energy for the different interaction types. 
The response matrices are shown graphically in Fig. H2] and numerically in 



App. Appendix A 



6. Conclusions 

Through a combination of fiducial, track quality, u u CC selection and 
kinematical cuts, an analysis has been applied demonstrating the power of 
MIND to detect and identify CC DIS interactions in the presence of real- 
istic reconstruction of the primary muon. The efficiency threshold currently 
lies between 3 and 4 GeV, and an efficiency plateau of 70% is reached at 
about 6 GeV. While improved sensitivity could be achieved by lowering the 
threshold, as mentioned in Sec. El this region would be dominated by quasi- 
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Figure 12: Response matrix in true/reconstructed neutrino energy for signal and back- 
grounds: i) signal efficiency, ii) backgrounds from CC, Hi) NC and iv) v e CC. 
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elastic and resonance interactions which have not yet been considered. This 
type of interaction should contain less hadronic activity and thus the low en- 
ergy pattern recognition efficiency should improve. An increased sample of 
successfully reconstructed events in this region should increase the efficiency 
of the golden channel analysis. 

The powerful rejection afforded by the inclusion of the hadronic energy 
and direction vector highlights the importance of good hadronic reconstruc- 
tion. While the hadronic energy has been reconstructed well in other similar 
experiments the direction vector requires careful consideration of both tech- 
nology and analysis to achieve the required resolution. In a future publication 
we will include a re-optimisied MIND design, within a GEANT4 framework, 
where we will also take into account low energy quasi-elastic and resonance 
interactions. 

Compared to the baseline MIND presented in |9_|, where perfect pattern 
recognition was assumed, these new results show some improvement. The 
aforementioned study considered the charge and NC backgrounds and in both 
cases the results presented here are of similar level. The signal efficiency curve 
reaches a plateau at 70% in the bin of 6 — 7 GeV. The corresponding curve 
in the previous study reaches approximately the same level at a similar or 
slightly higher energy depending on the particular analysis. Using this ef- 
ficiency curve, Neutrino Factory sensitivity studies were carried out in the 
context of the International Scoping Study (ISS) for a future neutrino facil- 



ity [29], demonstrating that a NF with two 50 ktonne MIND detectors at 
two different baselines has the largest #13 — 5c p coverage out of all possible 
facilities. This paper reinforces the conclusions from that study by showing 
that the pattern recognition and reconstruction of events in MIND do not 
degrade the selection efficiency for the oscillated signal. 
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Appendix A. Numeric summary of Response Matrices 

This appendix summarises the response matrices of signal and all back- 
grounds in bins relevant to an oscillation analysis. In all tables columns 
represent the true neutrino energy in GeV and rows the reconstructed en- 
ergy, also in GeV. The overflow bin in reconstructed energy represents all 
events with a reconstructed energy greater than the known maximum. 
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Table A. 4: Signal Efficiency response matrix; All values xlO 
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Table A. 5: CC background response matrix; All values xlO 
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Table A. 6: Neutral current background response matrix; All values xlO 
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Table A. 7: v e CC background response matrix; All values xlO 
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